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H ' In this work we propose a car cellular automaton model that reproduces the experimen- 

O tal behavior of traffic fiows in Bogota. Our model includes three elements: hysteresis 

^ , between the acceleration and brake gaps, a delay time in the acceleration, and an instan- 

taneous brake. The parameters of our model were obtained from direct measurements 
inside a car on motorways in Bogota. Next, we simulated with this model the flux-density 
fundamental diagram for a single-lane traffic road and compared it with experimental 
data. Our simulations are in very good agreement with the experimental measurements, 
C^ • not just in the shape of the fundamental diagram, but also in the numerical values for 

/— ^ ' both the road capacity and the density of maximal flux. Our model reproduces, too, the 

qualitative behavior of shock waves. In addition, our work identifies the periodic bound- 
ary conditions as the source of false peaks in the fundamental diagram, when short roads 
are simulated that have been also found in previous works. The phase transition between 
free and congested traffic is also investigated by computing both the relaxation time and 
the order parameter. Our work shows how different the traffic behavior from one city to 
another can be, and how important is to determine the model parameters for each city. 



Keywords: Cellular automaton models. Fundamental diagrams. Jams, Shock waves, 
jamming transitions. 



PJ ' 1. Introduction 

O ■ 

fi , Some years ago, Bogota was a city with heavy traffic congestion and a chaotic 

transportation system, just because it has 7 milhon inhabitants with more than 
55,000 taxis, 18,000 buses of different kinds, and a miUion of private cars roaming the 
rS I streets. Recent city administrations have tried to solve this problem by introducing 

j^ ' transportation strategies such as: a mass transportation system {Transmilenio) , 

almost 250 kilometres of bike paths, pedestrian bridges everywhere and restrictions 
on the use of private cars at rush hours {Pico y placa). Thanks these efforts, Bogota 
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has reduced the mean travel time in a 40%, accidents in an 80% and pollution in a 
50%. However, there are still many improvements to be done in the future. 

During the last 10 years, cellular automata models (CA) have been applied with 
success to trafSc simulations. These models are able to reproduce the macroscopic 
properties of highway traffic from the microscopic behavior of each car. The first 
model (STCA) ^'^, proposed by Nagel and Schrekenberg in 1992, is just based in 
the distance to the next car {gap=Ax) and the maximal speed, but leads to a 
quite realistic flow-density relation (fundamental diagram) and reproduces well the 
spontaneous jam and shock waves formations in highways ^''^. Later developments 
add other elements to STCA, like improved gaps including the speed difference 
to the car ahead ^'^, the speed at the previous time step **'^, and many other 
parameters. Some theoretical and practical studies have extended these models to 
two- and three-lane highways ^. 

The drivers' driving is very different from city to city, and a realistic traffic model 
should keep in mind the particularities of each place. Simulations and studies of 
this type have been carried out in cities like Portland, Los Angeles, Tokyo '^'^'S. 
However, there are not such studies performed in Colombia, even in Bogota. 

In this work we propose a car cellular automaton model that reproduces the 
experimental behavior of traffic flows in Bogota. Our model includes three elements. 
The flrst one are the gaps the driver uses to decide to brake (brake gap gapbrake) 
or accelerate (acceleration gap gapaccei)- They are, in general, different (hysteresis) 
and both depend on the speed. The second element is the time it takes the car to 
reach the next discrete speed value (retarded acceleration, tup). The last one is an 
instantaneous brake reaction (that) we have observed when the car ahead brakes. 
The parameters of our model were obtained from direct measurements inside a 
car that was running on Bogota's highways. With these parameters, simulations 
were performed to construct the flow-density fundamental diagram. This result was 
compared with experimental measures from Bogota's highways. The model was also 
used to compute shock waves, both in the free and congested regimes. Finally, we 
studied the phase transition between free and congested traffic (jamming transition) 
by computing both the relaxation time and the order parameter according to their 
definition in 14,15 'pj^jg }a,st study was performed on both the deterministic model 
proposed above and on the same model plus a probabilistic spontaneous-brake rule, 
as in STCA. 

The paper proceeds as follows. In section 2, we make a detailed description of 
our model, the values of the measured parameters and the vehicle rules. Section 
3 we show the simulation results and compare them with experimental data. Sec- 
tion 4 includes our study on the jamming transition. Section 5 contains the main 
conclusions and discussions of our work. Finally, appendix A and B describe in 
detail the experimental methods we used to obtain the model parameters and the 
experimental flux-density diagrams. 
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2. Our Model Description 



In our model, the highway is represented by a one-dimensional array of length 
L with periodic boundary conditions. Each site of the array is a cell of length 
2.5m, that is a finer discretization than the used in STCA model. Vehicles can 
only have integer velocity values, u=0, 1, ..., Wmax- We used frnax=7 and a speed 
unity Vunity — ^0^- This corresponds to time steps of tstep—0.9s, that is near to 
the usual value to driver's reaction time. These discretizations are usual for any 
CA traffic model, with the only difference that a vehicle occupies two consecutive 
cells, the length of a car (4.5m) plus the distance between cars in a jam (Im). 
Then, the maximal number of vehicles in the highway is given by N=^. At time t 
the n- vehicle is completely defined by: its position a;„(i), its velocity Vn(t) and its 
brake-light status, bn(t), which is 6„=1(0) when the driver brakes (or not) at the 
previous time step (t— 1), like ^^. The effective gap is defined as gap=Ax{t) + Av{t) , 
where Ax{t)=Xn+i{t) — a;„(t) — 1 is the number of cells empty to the vehicle ahead 
and Av{t)=Vn+i{t) — u„(t) is the speed difference to the car ahead. 



As already mentioned, our model includes three elements: the hysteresis between 
brake and accelerate gaps, the retarded acceleration and the instantaneous break. 
The three parameters gapbrakej gaPaccei and tup, are function of speed and they 
represent on the whole the drivers' driving. These parameters were experimentally 
found from inside a car on Bogota's highways (see appendix A) and are summarized 
in table 1. 



Drivers' driving in Bogota 
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Table 1. Drivers' driving from Bogota. 



Summarizing, all cars execute in parallel the following set of rules: 
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Rules 

• Compute its gap(t). 

• Read its parameters gapaccei, gciPbrake and t^p from table 1. 

- normal brake: If gap{t)<gaphrake, speed down to the maximal speed 
v{t+l) such that gap'^^^^^<gap{t)<gap'^^^^i, where (jaKrafce and gap'^^^^^ 
are the parameters at speed v{t + 1). In addition, let delay=0 and turn 
on brake lights (&„(i + 1)=1)- 

- If gap>gapacceh then 

* instantaneous brake: If gap(t)<gapaccei + 2 and the brake lights of 
the car ahead are on (bn+i{t)=l), let v{t+l)=v{t) — 1 (brake), turn 
on brake lights (6„(f + 1)=1) and let delay— . 

* accelerate: Else, turn off brake lights (6„(t + 1)=0) and 

• li delay==tup, let v{t + l)—v{t) + l (accelerate) and let delay=0. 
■ Else, let delay=delay + 1 and preserve v{t + l)=v{t). 

- Otherwise, let delay— 0, turn off brake lights (6„(t + 1)=0) and preserve 

v{t + l)=v{t). 

• Finally, move v cells ahead, 

x{t + 1) ^ x{t) + v{t + 1) . (1) 

The counter delay defines if t„p has been completed. The variable &„+i(t) defines 
the brake light status of the car ahead. The instantaneous brake rule represents the 
braking reaction we have observed when the car ahead also brakes. This reaction is 
observed for all distances but is just included in the gaps when gap < gapaccei- Thus, 
we have included it as an additional rule only it gapaccei < gcLpit) l£ gciPaccei + 2 
through a brake light on each car. 

3. Results 

The system starts with a initial configuration of N cars, with random distri- 
butions of speeds and positions and &j(0)=0 for all j. In order to prevent traffic 
accidents, as a previous step, we limit the speed values to the headway (wi„i<Aa;). 
Starting from this initial configuration, we measure the average velocity v{t) over 
all cars at each time step, t. After many time steps (t^oo), when the system 
reaches a stationary velocity state v{t)=v{oo), the flow is computed by the relation 
q=v{oo) ■ p=v{oo) ■ ^. The whole process is repeated 100 times for each density 
value, just to make statistics, and, so, the fundamental diagram is obtained. 

Our model is able to reproduces the phases observed in real traffic: free-fiow, 
synchronized, and stop-and-go. The shock waves and jams formations for these flow 
regimes can be observed in space-time plots (figure 1). 
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a) b) 

Fig. 1. Space-time plot for a road of L=1000: a) p=0.16 (p < pc), b) p=0.6 (p > pc)- Cars move 
from left to right and time runs from top to bottom. One can observe the schock waves and jams 
formations. If 7'ho < p{qmiix){rho > p((?max)), these formations are moving forward (backward) 



Before comparing the simulation results with the experimental data, we study in 
detail the finite size effects. For small systems, appears in the fundamental diagram 
a false peak of maximal flow, which disappears when the system is sufficiently 
large. Figure 2 shows the fundamental diagram for different system sizes. One 
can observe that the false peak appears around p—0.12 (in the synchronized traffic 
phase), and it has completely disappeared for L=2000. These false peaks have 
also been observed in many other models i^'^i. In ^^ they are interpreted as the 
coexistence of two phases of traffic flow in a dynamic equilibrium. In contrast, we 
have found that these false peaks are due to extraordinary configurations that exists 
only in small systems. Due to the periodic-boundary conditions, it is possible to 
think a configuration where all cars have the maximal speed and all have the same 
gap, equals to that maximal speed. 

In these configurations the system does not relax, but remains forever, with 
a mean velocity that is larger than the average velocity for the relaxed system. 
Thus, they push up the values of (^(cxo)), where the () denotes the average on many 
realizations, generating the false peak. In conclusion, the peak should be located 
at densities p=I/i'niax=0.14, in good agreement with the computational results. 

Since our model takes into account three elements, we were interested in looking 
at the effect of including each one of them, one by one, namely: the hysteresis be- 
tween brake and acceleration gaps, the retarded acceleration and the instantaneous 
brake. 

Figure 3 shows the fundamental diagrams obtained from our model with (a) 
brake and acceleration gaps only, (b) gaps plus retarded acceleration and (c) all three 
elements. They are compared with (d) the fundamental diagram from the STCA 
model with our time and space discretization. One notes small differences between 
the fundamental diagrams corresponding to STCA and the most elementary model. 
This means that the hysteresis doesn't have any relevant effect on the fundamental 
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Fig. 2. Eflfccts of small system sizes (L) into the simulated fundamental diagram. One observes 
the false peak of maximal flux due to the finite boundary eonditions and how it disappears when 
L=2000 



diagram. Its consequences, however, should be seen in the behavior of jams and 
shock waves, but they will be not investigated here and they would be theme of 
future work. 

When the retarded acceleration is included, two effects emerge, namely: a de- 
crease of the maximum flow and a its shift towards larger densities. Finally, the 
instantaneous brake effect shows up slightly into the congestion region. 



STCA Model 

Our Model + Hysteresis 

Our Model + Retarded Acceleration 

Total Model 




Fig. 3. Effects that have each one of the elements of our model on the shape of the fundamental 
diagram. Also is shown the fundamental diagram of the STCA model with our time and velocity 
discretization. 
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To validate our model we performed a comparison with experimental measure- 
ments. Figure 4 compares the flow-density diagram from our model with mea- 
surements over Bogota's highways on broad density ranges, both as a) dispersed 
data and b) averaged data. A second alternative, consisting of a STCA model plus 
retarded acceleration , is also included. 

On one hand, one observes that both simulations are in good agreement with 
the experimental data, especially if we compare the numerical values, simulated and 
measured, of the road capacity and the density of maximal flux. 

P(9max)exp = 0.33(3) Qmaxexp = 1.43(6), (2) 

P{qma^)sim = 0.33(4) ^max sjm = 1.320(4). (3) 

The difference in the scatter between simulated and experimental data shows up 
(at least in part) because the first ones are averages over much more cells than the 
second ones. The two models have different behaviours just in the zone of high 
congestion, were our model predicts lower fluxes. This is the effect of hysteresis 
that is expected to play a role in jam formations. On the other hand, these values 
of road capacity are slightly larger than those measured in other countries ^ . This 
suggests that the traditional Bogota's aggressive driving makes the traffic fiux more 
efficient. The price is, however, one of the most high rates of fatal victims on the 
world ^^ (in fact, one out of six victims of violent causes in Colombia dies in a car 
accident "'^^). For a recent work including agressive drivers see ^^. 

Cellular automata models for traffic flow exhibits sometimes a phase transition 
from a free-flow phase to a congested phase I4,i5,i6 -^u most cases this transition is 
flrst order, but some models, like STCA, shows a second-order phase transition at 
a single point in the phase space ^'^ . The relaxation time and the order parameter 
are the quantities which characterize this transition. 

To compute the relaxation time we employed the definition of Csanyi and Kertesz 
^^ . As they do, we start from a configuration random initial positions and zero speed 
for all cars. Then, the relaxation time is computed as 

/■OO 

r= / [Tmn{v*{t),{v{^))}-{v(t))]dt , (4) 

Jo 

where v*{t)^t denotes the speed a car obtains at time t when there are no cars 

ahead (free acceleration). 

As we already know, a characteristic feature of a second order phase transition is 
the divergence of the relaxation time at the transition point. Figure 5 shows how the 
relaxation time has a maximum at a pseudo-critical density /?c=0.33, which is the 
same value of maximal-fiow density one can read from the flux-density fundamental 
diagram (figure 4), as expected for a deterministic model ^^. 

The typical order parameter for traffic cellular automata is ^^. 
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Fig. 4. Comparison of the flow-density diagram for tlic simulated models with the experimental 
measures. Both simulated models, our model and de STCA model plus retarded aceeleration, are 
in good agreement with the actual data. Above: dispersed data. Below: averaged data. 
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Fig. 5. Relaxation time r near the transition density pc for different system sizes. 



Since a car in our model occupies two consecutive cells, we redefine the order 
parameter as 



1 



(6) 



Like , figure 6 shows that the order parameter exhibits a continous transition. 
The situation is quite similar to the behavior of the order parameter in finite sys- 
tems. That is, the order parameter have small values for small density, and around 
the Pc it begins to have a non-zero value. 
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Finally, we wanted to explore what occurs when we add the randomization rule 
of the STCA model. This is, 

• Randomization: If after the above steps the velocity is larger than zero {v > 0) , 
then, with probability p, let v{t)=v{t) — 1. 

Figure 7 shows the fundamental diagrams for differents values of p. One observes 
that when the noise p increases, the flow decreases and the maximum towards 
smaller densities. In addition, for larger densities, the randomization rule causes an 
early collapse in the system. 




Fig. 7. Fundamental diagram for different values of p. This figure corresponds to our model plus 
randomization rule of the STCA model. 
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Fig. 8. Behavior of the order parameter for Our model with a randomization rule, a) Order 
parameter for different values of p and b) Scaling plot for the order parameter (excluding p=0). 
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Figure 8 shows the order parameter behavior for difFerents values of p. It con- 
verges smoothly to zero for densities smaller to pc- To the right hand, we performed 
the respective scaling. The following form summarize the simple scaling: 

M(p)-m(p+|) . (7) 

Herein A is the shift of the transition density compared to the same value for 
p=0.125, because it was impossible to make the scaling for the deterministic case. 

4. Conclusions 

In this work, we present one of the first cellular automaton models that applies 
to the reality of Bogota. Their objective is not just to propose a first model, but 
also, to investigate which models apply to our city. The work includes also measures 
of the driving parameters of a car in Bogota and of the fundamental diagram of a 
highway. 

Our model includes three elements, to summarize are the drivers' driving: hys- 
teresis between brake and acceleration gaps, retarded acceleration and instantaneous 
brake. 

The simulation results are in a good agreement with the experimental mea- 
surements. That is, the simulated fundamental diagram reproduces successfully 
the shape of the fundamental diagram and the numerical values for both the road 
capacity ((/max) and the density of maximal flux (p(9max))- The obtained values 
are: 



P{qmax)exp = 0.33(3) ^max exp = 1.43(6) (8) 

P{qma^)sim = 0.33(4) ^max sim = 1.320(4). (9) 

For small sisteni sizes we have found that the false peak of maximal flux in the 
fundamental diagram is due to some extraordinary configurations. In these config- 
urations the system remains in a not-relaxation state, with a mean velocity larger 
than (u(cxd)). Finally, these extraordinary configurations only exist for systems with 
L < 2000. 

The hysteresis element in our model doesn't have any relevant effect on the 
fundamental diagram. If we compares it with the STCA model, the differences 
appears only into the congestion region. However, we believe that the spontaneous 
formations of jams and shock waves should be affected by this elements, and this 
is an interesting area of future work. In contrast, the retarded acceleration is much 
more important. It fixes the real shape of the fundamental diagram and the char- 
acteristic values of gmax and p((7max) • Moreover, an STCA model plus the retarded 
acceleration is enough to reproduce the fundamental diagram we found in Bogota's 
highways. 

By looking at the phase-transition, we found qualitatively the same behavior 
of both relaxation time and order parameter as shown in I4,i5 q^-^^^ observes a 
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maxiniuin of the relaxation time t near the density transition pc, specially for the 
deterministic case Pc—piqmax)- The divergence is more marked when the system 
size increases. The behavior of the order parameter shows a smooth convergence to 
zero for small densities. 

Summarizing, both our model an the STCA model plus retarded acceleration 
reproduces the fundamental diagram for Bogota. They coul be used to perform 
more complicated simulations in future works, like semaphorized, intersections and 
even the Bogota's net of main highways. It also remains to study the jam formation 
in both models. 
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Appendix A : Our Model Parameters 

By installing a videocamera inside a car, we taped the driving of a driver on 
the highways in Bogota. The videocamera was calibrated with the plates of the 
other cars, therefore we made a chart, which relates the longitude of the plates with 
the distance to the car in front. Ax. During the trips, the copilot registers the car 
speed the velocity and if the driver brakes, accelerates or preserves this speed. With 
the films on hand, one can know at any time the vehicle speed, its state (braking, 
acceleratimg or preserving) and the distance to the car in front (Ax). The relative 
velocity Av is calculated as 

Axa - Aa::i 
Aw = , (1) 

where Axi is the usual Ax, and Aa;2 is the distance to the car in front measure 
after two seconds. 

To include the retarded acceleration it is neccesary to measure the parameter 
tup- For this purpose, the car was accelerated from 0^ to 100-^ and we measure 
the time, it takes to reach the next discrete speed value, in steps of (10-^). 

Appendix B : Experimental Measurements 

The experimental data were obtained by capturating on tape (from a pedestrian 
bridge), the traffic flow in a highway segment. The 30*'* avenue is perhaps the most 
important avenue of Bogota, because it communicates directly the south with the 
north of the city. At daily rush hours the 30*'* avenue is not able to cope with the 
demand and bored traffic jams are generated. We chose for our measurements the 
south-north high-speed lane of the 30*'' avenue between 53*'* street and the Campin 
Football Stadium. 

This is a two-lane sector and have a length of L =169(5)to«62(2) cells, i.e, a 
maximal number of cars iVmax=31(l). 
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Fig. A.l. Digitized pictures of the tape when the parameters of the model were measured. 

We carried out the following process: 

• At each time, the density of the system is computed by p= jy i where N(t) 
is the number of cars over the highway sector. 

• The velocity of the system is the average velocity over the N cars. It is cal- 
culted as an aritmethic mean, v{t)—jj J2j=i ^j(^)i where Vj{t) is the average 
velocity of each car over the highway sector. 

• Finally the flow is calculated as q~p ■ v{t). 

With this process we obtained the experimental data to perform the flux-density 
fundamental diagram. 




Fig. B.l. Digitized picture of the tape, when the experimental fundamental diagram was measured. 
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